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1 Introduction 

This article presents a tutorial introduction to a recently developed real-time 
renormalization group method [Q . It describes nonequilibrium properties of dis- 
crete quantum systems coupled linearly to an environment. We illustrate the 
technique by a simple and exactly solvable model: A quantum dot consisting of 
a single non-degenerate level coupled to two reservoirs. The article is intended for 
advanced students. Besides elementary quantum mechanics and statistical me- 
chanics, it requires knowledge of second quantization and Wick's theorem. The 
latter topics can be learned easily from standard textbooks, see e.g. Ref. 

Renormalization group (RG) methods are standard tools to describe various 
aspects of condensed matter problems beyond perturbation theory . Many im- 
purity problems have been treated by numerical RG with excellent results both 
for thermodynamic quantities and spectral densities These RG techniques, 
however, cannot describe nonequilibrium properties like the nonlinear conduc- 
tance, the nonequilibrium stationary state, or the full time development of an 
initially out-of-equilibrium state. To address these aspects we present here a per- 
turbative RG method, formulated for strongly correlated quantum systems with 
a finite number of states coupled linearly to external heat or particle reservoirs. 
Examples are: spin boson models, molecules interacting with electrodynamic 
fields, generalized Anderson-impurity models, quantum dot devices, magnetic 
nanoparticles interacting with phonons, etc.. Fundamentally new, we generate 
non-Hamiltonian dynamics during RG, which captures the physics of finite life 
times and dissipation. Furthermore, no initial or final cutoff in energy or time 
space is needed, i.e., large and small energy scales are accounted for correctly 
like in fiow-equation methods Although correlation functions can also be 
studied, physical quantities like spin and charge susceptibilities or the current 
can be calculated directly without the need of nonequilibrium Green's functions. 

The purpose of our RG technique is to describe quantum fluctuations which 
are induced by strong coupling between a small quantum system and an envi- 
ronment. There are several recent experiments which show the importance of 
quantum fluctuations in metallic single-electron transistors Q and semiconduc- 
tor quantum dots |8| (see |^ for an overview over theoretical papers). Due to 
the renormalization of resistance and local energy excitations, anomalous line 
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shapes of the conductance have been observed, which can not be explained by 
golden-rule theories. For applications of the real-time RG to these cases we refer 
to Rfs. 

Here we want to illustrate the method by an exactly sovable model, namely 
a quantum dot with one non-degenerate state with energy e coupled to two 
reservoirs (r = L,R). The Hamiltonian H = Hn + Hq + Ht consists of three 
parts, corresponding to the reservoirs, the dot, and tunneling 

^^fl = ^ -ffr ^ l^kra\r°-kr , (1) 
r kr 

Ho ^ec^c, (2) 

Ht=J2 {^ralc + T^c^akr} ■ (3) 

kr 

Nonequilibrium is taken into account by describing the electrons in the reser- 
voirs by Fermi distribution functions with different electrochemical potentials 
(Ur- Tunneling is switched on suddenly at the initial time to, i.e. initially the 
density matrix p{to) = po decouples into an equilibrium part for each reservoir, 
Pr = Z~^e~^'^^'^~^'^^^\ and an arbitrary initial distribution p{to) — Po for the 
dot 

Po = PoPres = PoPLPres ■ (4) 

The aim is to calculate the time evolution of the reduced density matrix of the 
dot, p{t) — Trresp{t), and the tunneling current {Ir){t) flowing from reservoir r 
to the dot. The tunneling current operator is given by |^ 

Ir = -eNr = {^rO-irC - T* Ukr^ . (5) 

The solution of the above quadratic Hamiltonian is trivial since all degrees 
of freedoms can easily be integrated out. Doing this by using Wick's theorem 
for all field operators within the Keldysh formalism, one can easily solve the 
full nonequilibrium problem However, except for having solved a special 
and almost trivial problem, we would not have gained anything for solving more 
general problems of dissipative quantum mechanics. Usually, the local system 
can not be integrated out due to interaction terms or spin degrees of freedom. 
Therefore, we try to proceed differently. We will only integrate out the reservoirs 
and keep the dot degrees of freedom explicitly. This is always possible for an ef- 
fectively noninteracting bath. As a result, we get an effective theory in terms of 
the local degrees of freedom, expressed by a formally exact kinetic equation for 
the reduced density matrix of the dot. For the special Hamiltonian we 
solve this equation exactly. Furthermore, we will also develop a renormalization 
group method to solve the kinetic equation. We show that the RG equations 
describe the same exact solution. The important point is that both steps, i.e. 
(a) setting up the kinetic equation, and (b) setting up the RG equations, are 

^ Throughout this work we set h = ks ~ 1 and use e < 
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not specific to tlie above Hamiltonian but can be applied to any discrete quan- 
tum system coupled linearly to an environment. The only difference is that, for 
most problems, the resulting RG equations have to be solved numerically. In 
conclusion, the above Hamiltonian serves as a test example to illustrate the RG 
technique and to demonstrate that it is well-defined and useful. 

2 Diagrammatic language 

2.1 Diagrams on the Keldysh contour 

We start by introducing some convenient notations. The index fj, ~ rjr labels 
the possible tunneling processes between reservoirs and dot. 77 = ± indicates 
tunneling in/out, and r — L, R specifies the reservoir. We define the following 
reservoir and dot operators 

i-r = ^'■"fer ' ■^+^ ^ '^rO'kr , (6) 

k k 

, g+r = ■ (7) 

We denote the two possible states of the dot by s = 0, 1 with energies Eq ^ 
and El = e. The Hamiltonian (|l|)-(§) becomes 

Hji = Y likra\r°-kr , (8) 
r 

s 

Ht Qt^it. X! i^+rJ+r + j-rg-r} , (10) 

/J r 

where the symbol : . . . : denotes normal ordering of Fermi field operators but 
without sign change when two operators are interchanged. The current operator 
for the left reservoir is given by 

L = ^ : : - V = -ievgt^SrL , (a^ = Vr) , (11) 

and a corresponding equation for In. 

The time evolution of an arbitrary observable a follows from the von Neu- 
mann equation 

{a){t) = Trapit) = Tr a e-'"'^*-'°'> poe'"^'-'°^ (12) 
= Tre*^(*-*°) a e'^^^*"*" Vo , (13) 

where, in the last step, we have used cyclic invariance under the trace. To get 
a matrix element of the reduced density matrix of the dot, p{t)ss' = 
need a — \s'){s\. For the current, we take a ^ II = ■ 
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To integrate out the reservoirs, we expand the propagators in tunnchng and 
apply Wick's theorem to the reservoir degrees of freedom. We introduce the 
interaction picture 

b{t) = SHR+Ho)(t-to) ^-i(HR+Ho){t~ta) ^ ^^^^ 

and obtain 

mit t \ f \ - if* dt' HtU') , , -i V dt' HtH') 

g«_f/(t^to)^g-»//(t-io) ^j^g Jto ^^'a{t)Te •''0 (15) 

where T and T denote the time-ordering and anti-time-ordering operators, re- 
spectively. Inserting ( [l5| ) in (|l^) and expanding in Ht gives a series of terms 
which we visualize diagrammatically, see Fig. |ll The upper (lower) line corre- 
sponds to the forward (backward) propagator. The diagram shown corresponds 
to the following expression 

i'{-ifTrHT{h)HT{t2)a{t)HT{ti)HT{U)pQ . (16) 

We see that the operators are ordered along a closed time path (Keldysh con- 
tour), as shown in Fig. |l} 



— • > •- 



Po= Po Pn 



-• ^ 




iHT(t3) iH^Ct,) 

H \ \ 1 h 



to t4 t3 t2 ti t Time 

Fig. 1. Diagram corresponding to Eq. (p^). 



The next step is to insert Ht — '■ '■ for the tunneling Hamiltonian 
and Po — PoPres for the initial density matrix. We use the short-hand notation 
9i — 9fj.i(ti) and ji = j^i{ti), and decompose each diagram into a dot and 
reservoir part. If a = |s' >< s| is a dot operator, we get from ( |l6|) 

i'^i-if X! {Tro 9392 a{t) 9194 Po} {Trres hj2jljA Pres} , (17) 

whereas, if a = = : i^j^ : is the current operator, we get 

*^(-«)^ X! {TrQg?,92i^,{t)9l9'iPo}{Trreshj2jt,{t)jij4Pres}- (18) 
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Here, Ttq {Trres) denotes the trace over the dot (reservoir) degrees of freedom. 
The reader can convince himself very easily that this factorization does not imply 
any additional minus signs from commutation of Fermi operators. The reason is 
the quadratic form of the tunneling Hamiltonian. 

The trace over the reservoirs can be calculated easily by using Wick's theorem 
As a result, we can decompose any average over products of reservoir field 
operators into a sum over products of pair contractions. Denoting by (. . .) the 
average over the reservoirs, we get for the reservoir part of ([l7| ) 

I — I I — I I — Fn I I Fn I 

= (J3j2)(jlj4) - (J3jl)(i2j4) + {hj4){j2jl) ■ (19) 

Each pair contraction corresponds to an equilibrium average over two reservoir 
field operators. If two contractions intersect proper minus signs have to be taken 
into account due to Fermi statistics. A pair contraction {jp,{t)jfj,'{t')) depends on 
the relative time t — t' and is only non-zero for fi' = ft, with p, = —rjr (if /i = i^r). 
We define 

7p(i) = (j^(i)jM) . (20) 
and get, by using the definition (||) 

aL(i)afer) for rj = + 

v = - 



I (afcr(i)aL) for 



= ^ y dErr{E)e""''f'{E - fir) (21) 

c^£l f dEe'J^^-t'-^/^e''^^'f^{E) . (22) 
27r J 

Here we have defined 

Er = EriE) = 27rJ2\Trf^{E - ekr) = 2'K\TrfNriE) , (23) 

with Nr{E) being the density of states of reservoir r. We define /+ ^ f , = 
1 - /, and ff{E) = P{E - /i^), with f{E) = l/(exp(/3£;) + 1) being the Fermi 
function, and (3 = 1/T the inverse temperature. We take the density of states 
Nr{E) independent of energy and regularize the integral (|2^) by introducing 
an exponential convergence factor e''^^"'''')/^. Here, D corresponds to a high 
energy cutoff. In the end, we will send D — > oo. Performing the integral ( p2|) 
gives 

~ 2;5ih^h[^^(r^77i5)7/3] ■ ^^^^ 

Furthermore, we define 

int)^lM^ (25) 
and note the important property 

hm {j+{t)+j-it)}=rrS{t) , (26) 
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Fig. 2. Diagram after tlie application of Wick's theorem. On tlie propagators tiie states 
s = 0, 1 of tiie dot are sliown. At the vertices we have indicated the index /i = rjr. 



which follows directly from (p2[). 

We indicate the pair contractions diagrammatically by connecting the corre- 
sponding vertices by a dashed line, see Fig. ^. The remaining part for the dot 
degrees of freedom still remains, see Eq. (pT|). We calculate this part by inserting 
intermediate states of the dot between the operators. These states are indicated 
in Fig. ^ between the tunneling vertices. The reservoirs are already integrated 
out, so the tunneling vertices correspond to the dot operators g^,- (with an ad- 
ditional factor =Fi for a vertex on the upper (lower) propagator). Between the 
tunneling vertices we have the free time evolution of the dot, i.e. for a propaga- 
tion of state s from t2 to ti, we get an exponential factor e"*^'^*!"*^). 

Our diagrammatic language provides an effective description in terms of the 
dot degrees of freedom. The presence of the reservoirs is reflected by the retarded 
coupling of the tunneling vertices by the free Green's function of the reservoirs. In 
particular, this means that the forward and backward propagator are n o lo nger 
independent but are coupled by reservoir lines. We will see in section 3.1 that 
this leads to rates in a kinetic equation for p{t). 



2.2 Superoperator notation 

In this section we will replace the double-propagator diagrams on the Keldysh 
contour by a convenient matrix notation. This provides a very compact and 
analytic way to express diagrams by formulas. 



Po • • • • • • • 4 • 10x01 

00 10 00 01 00 01 11 01 00 

Fig. 3. Diagram in the double-state representation. It results from Fig. |^ by taking the 
upper and lower line together to one single line. 



Instead of considering two propagators and specifying the dot states on each 
propagator separately, we can formally take both propagators together to one 
line and specify the state on this new propagator by a double dot state (s, s'), 
see Fig. ^. Here, the first (second) state corresponds to the state on the upper 
(lower) propagator. By this trick we have lost the information wether a tunneling 
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vertex lies on the upper or lower propagator. To recover this we add to the index 
II of the tunneling vertex operator an additional index p = ±, which indicates 
wether acts on the upper (lower) propagator. The new vertex is denoted by 



Gf, with matrix elements 



(27) 
(28) 



We have included the factors =Fi for a vertex on the upper (lower) propagator 
into the definition of G^. The free time evolution between the vertices is given 
by a factor e"*(^=~^s')^*i~*='', where (s, s') indicates the double state on the line. 



This can be written in operator form as (e *^o(*i with 



-iLo(*i-t2)'| 
(io)sis'j,S2s!, — '5siS2'5s'jS^(-E'si — Eg 



(29) 



Finally, a contraction connecting a vertex G^, at time t' with a vertex G^ at 

time t > t' , is denoted by 7^^'(^ — t')- Using the definition ( po|) and the fact 
that operators on the lower propagator act always later than those on the upper 
propagator, we obtain 



(JpWjp)=7a.W for p' = + 



2/3sinh[7r(p't-V^)//3] 



(30) 
(31) 



where we used the result ( p^ in the second step. 

The double-state matrices Lq and Gj^ are called superoperators in the sense 
that they act on single-state matrices, i.e. on ordinary operators. If b is an 
ordinary operator, we can define Lq and G^^ by 



Lob = [Ho,b] , G+& = -i9f,b , G^ 6 = 165^ . 



(32) 



Within the superoperator notation and using a = \s'){s\, an arbitrary dia- 
gram for pss' (t) can be written as 



(33) 




The time dependence of Gi — G^'.{ti) is defined by 

GP(t) = e»^«(*-*«)GPe-*^°(*-*'') 



(34) 



All operators Gi . . . G„ are coupled in all possible ways by reservoir pair con- 
tractions, as indicated in (p3|). Implicitly we assume summation over . . 
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and pi . . .pn, together with the integration over the time variables ii . . . t„ with 
t > h > t2 > . . . > tn > to. 

If a = /l = : ifj.j^ '■ corresponds to the current operator, we get 



In comparism to (|33|), we need an additional pair contraction to the current 
vertex. In order to treat the boundary vertex as well within the superoperator 
notation, we define a superoperator by 

I+b = t^b/2 , (36) 

with matrix elements 

Using cyclic invariance under the trace, we get for (pq) 




(38) 

where /( — Ifi_(t), and the interaction picture is defined by 

IP{t) = IPe-'^''^'-*''^ (39) 



Eqs. ( |33| ) and (|3^) are the central result of this section. They relate the 
reduced density matrix of the dot and the average current to diagrammatic 
expressions in a very compact and analytic way. It turns out that the usage of 
superoperators simplifies the notation considerably. We will see in sections ^ and 
^ that the derivation of kinetic equations and renormalization group equations 
is very transparent in this language. However, one should always keep in mind 
that the usage of superoperators is only a formal trick to find a convenient 
matrix notation. Therefore, we have set up the diagrammatic representation in 
terms of the Keldysh contour first, and then, in a second step, introduced the 
superoperators. Of course it is also possible to start directly with superoperators 
[Ej, which provides a more compact and shorter way to arrive at Eqs. (|3^ ) and 
(p8|). However, for pedagogical and physical reasons, we did not proceed in this 
way here. The diagrams on the Keldysh contour reveal better that there are 
different kinds of terms which have to be distinguished very carefully from a 
physical point of view. Reservoir lines connecting the upper and lower propagator 
correspond to rates, they change the state of the dot simultaneously on the upper 
and the lower propagator. This describes a transition from one diagonal matrix 
element of the reduced density matrix p{t) to another one. Such processes can 
not be expressed on a Hamiltonian level and lead basically to the physics of 
dissipation. In contrast, reservoir contractions which connect vertices within the 
upper or lower propagator describe renormalization and broadening of levels. 
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3 Kinetic equation 
3.1 General approach 

In this section we will derive a self-consistent equation for the reduced density 
matrix p(t) of the dot, together with an expression for the average current. To 
achieve this it is essential to distinguish in Eqs. and (Q) between connected 
and disconnected parts. From (|33|) we see that any diagram for p(t) can be 
written in the form 



_^ g-»Lo(t2„-2-t2„-l) (A2„_lG 



-iLo(t-2-t3) 



(A3G...GB4) con • • ■ 

.GB2„)con e-^^^^'^-^o^Po ■ (40) 

Here, {AiG . . . GBi^i)con denotes a sequence of vertices between ti^i and ti 
which are coupled by pair contractions in such a way that any vertical cut 
between ti+i and ti will cross some contraction. We define such a block as a 
connected part of a diagram. E.g. Fig. ^ shows a sequence of three connected 
blocks. The boundary vertices A and B arc identical to G, i.e. = B^ — Gf^, 
but the interaction picture is defined differently 



(41) 



The reason is that we want the connected part {AiG . . . GBi+i)con of a diagram 
to depend only on the relative time argument ti — t^+i. Furthermore, we distin- 
guish the boundary vertices and B^ from G^ since, within the renormaliza- 
tion group procedure developed in section ^, the boundary vertices renormalize 
differently. 

We define the sum over all connected diagrams between t' and t by the kernel 

S{t - t') 



S{t - t') 



{AtGiG2 ■ ■ ■ G2nBt 



(42) 



We note the important property that the kernel S{t — t') is independent of the 
initial time to since all exponential factors e^'^"*" arising from the interaction 
picture cancel within S. Thus, in order to calculate we can set io = in the 
definition of the interaction picture of A, B, and G, see Eqs. ( p^ and (|4l|). 
Using the definition (^ ) for 2J in (^0|), we obtain 



Pit) 



E 



dti / dt2 ■ 

to J to 



dt2 



-'^o(t-ti) jj^^^ _ ^-^Lo(t2-ts) _ 



-iLo{t2 



(43) 



Differentiating with respect to time gives the kinetic equation 



Pit) + iLop{t) - /,*^ dt'S{t - t')p{t') 



(44) 
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Since the r.h.s. of this equation is a convolution in time space, we can formally 
solve this equation in Laplace space. We define the Laplace transform by 

/•oo 1*00 

p{z)= dte''*p{t) , S{z)^ dte'''U{t) , (45) 

J to Jo 

and get from (Q) the solution 

Z ^ Lq — tZj[Z) 

The time dependence p{t) follows by reversing the Laplace transform 

00 

pit) = lim — / dz e-"'p{z) = lim — / duj e-'"*p(w + ir]) . (47) 



We remark that p(z), defined by (45), is analytic in the upper half plane since 
p{t) will approach a stationary value for t —>■ 00. Thus, within the integration 
region of ( |4^ ) the integrand is well-defined. We see that the integral kernel S{z) 
is the central object which has to be calculated. The full time evolution of p{t) 
out of an arbitrary nonequilibrium state can be obtained once S{z) is known for 
all z — cj + iO+. The calculation of S{z) will be the subject of the renormalization 
group approach described in section 0. 
The stationary state is defined by 

Pst = lim p{t) — —i lim zp{z) . (48) 

Multiplying ( p6| ) by z[z — Lq ~ iS{z)] and taking the limit z iO+, we see that 
the stationary state is the eigenvector of Lq + iIJ{iO'^) with eigenvalue zero 

[Lo + S{iO+)]Pst ^ ■ (49) 

The density matrix p(t) is hermitian. In Laplace space this is equivalent to 

p{z)l, = p{-z*)s's ■ (50) 

This follows from the solution ( ^6| ) by using the symmetry relations 



where p — —p. They follow directly from (|27|), (|28D, (p9|), and the hermiticity of 
the Hamiltonian. The consequence for the kernel ( |4^ ) is 

, • (52) 

Using these relations together with the hermiticity of the initial density matrix 
Pq, we find directly (E^) from (E6|). 
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To prove conservation of probability Trop{t) — 1, we first note that 

^(Lo)..,. = = Y.(Bl)ss.. = Y.(Gl)ss,.. = . (53) 

s sp sp sp 

Applying this to (|2|) together with the fact that the contraction connected to 
the boundary operator does not depend on p, we find the same property for 
the kernel 

s s 

Applying Ttq to the kinetic equation (^), and using the properties (^3|) and 
(p4[), we find d/dtTrop{t) = which proves conservation of probability. 

To calculate the current we proceed analogously. From (|3^) and (J^j we see 
that any diagram for the average current can be written as 

TrQ{ItG...GBt^),onPit') . (55) 

The first connected block contains the current vertex. The sum over all connected 
diagrams containing the current vertex is denoted by Si{t — t') 



{ItGlG2 ■ ■ ■ G2nBt 



(56) 



The only difference to (^) is that the boundary vertex A^^ has been replaced by 
the current vertex /p. From (^5|) and ( p6| ) we find 



{iLm - E;.(: VJm :) = Si dt'Tro Si{t - t')p{t') 



(57) 



In Laplace space we get 

{lL){z)=TroSi{z)p{z) , 
and the stationary solution follows from 

{lL)st = Tro Ei{iO'^)pst ■ 



(PY , 



- (TP 



(^/i)s'j^si, 



and 



(58) 



(59) 



The expectation value of the current is real, i.e. {Il){z)* = {Il){—z*). This 
can be seen from the solution (58) by using the symmetry relations 



(60) 



(61) 
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3.2 Exact solution 

For the model of a single non-degenerate dot state, the kinetic equation and 
the current formula can be solved exactly. We derive this solution here and will 
check in section ^ that the renormalization group equations describe the same 
solution. If the reader is not interested in technical details, he can find the final 
results in Eqs. ( |65| ) and (|7|), and can proceed to the next section. 

There are only two possible dot states s = 0, 1. We denote by s the conjugate 
state, s = 1(0) if s = 0(1). Furthermore we define p — —p and p, — —rjr if /x = rjr. 
The following three properties are needed for the following 

J2p{I^ss,..^0 (62) 

ps 

(63) 



(64) 

The time arguments of the interaction picture are not written explicitly. The 
proof can be found in the appendix. 

For the reduced density matrix p{t) and the current (/i)(t), we need the 
kernels ^(2) and Si{z), see Eqs. ( ^6| ) and (|58|). Due to particle number conser- 
vation of the total Hamiltonian, the reduced density matrix p(t) stays diagonal 
if the initial density matrix po is diagonal. We therefore need only the diagonal 
matrix elements S{z)ss,s's' and Z'/(2;)ss,s's' (the nondiagonal elements can also 
be calculated but do not contribute to the current due to Tvq in Eq. (^)). 

The kernels are defined in ( ^2| ) and (|5^). Using property ( |62| ) we find 



^l{z)ss,ss — ^l{z)ss,ss ■ 



As a consequence, we only have to calculate the matrix element (11,00) of 
the kernels (the matrix element (00, 11) is analog and we quote the result at the 
end). 

We call an intermediate double-propagator state ss "even", and a state ss 
"odd". The vertices G and / change the parity. We want the matrix element 
(11,00), i.e. we start with an even state 11 from the left. Thus, the state after 
the first vertex of the kernels in Eqs. (^ ) and (|56| ) is odd, and we can apply 
properties (|6^ ) and ( |6^ ) for the pair of the second and third vertex. We see that 
the only possibility is that these two vertices are connected by a pair contraction. 
After these two vertices we are again in an odd state and can proceed in the 
same way. Finally we find that the only nonvanishing diagrams for S{z) are 

AGG GG ... GGB , (66) 



ps I 

E(g^gJ),,,. = o 



pp' 



lim YiGlG^'] 



(65) 
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and analog for Sj{z) by replacing A ^ I. 

I — I 1 — 1 

We denote the sum over all sequences of GG -blocks in (66) by 11 {t), where 
t is the time difference between the boundary vertices A and B. n{t) can be 
calculated analogously to section 3.1 where we resummed sequences of iT-blocks 
with the result (^6|) for the reduced density matrix in Laplace space. Thus we 
obtain 



n{z) 



with 

a{t) 

The kernels follow from 



GG 



z — Lq — ia{z) ' 
^iZWGle-^'^^'Gl, . 



(67) 
(68) 

(69) 
(70) 



Using all definitions together with (|2^) and the fact that Lq is a diagonal matrix, 
wc find in the limit D oo 



m -iZ'^t)Gin{t)Gl, , 



o'(i)io,io 
o'(z)io40 

n{z)io.iQ 



n{t) 



10,10 



-r/2 , 

i 

z - (Lo)i040 - ia{z) 

-iet -rt/2 

e e ' , 



-rs{t) , 



10,10 



■ir/2 



(71) 



with "f^ = 7r ^'id r = X^r-^f Furthermore, using the symmetry relations 
(|l]) we get iI(i)oi,oi = ^(Oio,io- Using these results together with 77(t)* = 
Jr{—t), we can evaluate (|69|) and ([70|), and find 



^(t)ii,oo = j+{t)n{t) 10.10 + h.c. 



(72) 



i:/(i)ii,oo = (e/2)7+(t)7T(i)io,io + h.c. = (e/2)7+(i)e-^^*e~"/2 ^ ^.0.(73) 



Using (E2|) and (p5|) we finally get in Laplace space 



^(^)ii,oo= ^V/^'^E.^r/.l^^)! 
^/(^)ii,oo = ||/di?rL/i(i?){ 



E~e+z+ir/2 
1 



E~e-z-ir/2 
1 



E-t+z+ir/2 E-t-z-ir/2 



(74) 



By an analog calculation one obtains the matrix element (00,11) by replacing 
fr^^~ fr and changing the sign for the current kernel. 
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4 Renormalization group 

In this section we will develop a renormalization group technique to calculate the 
kernels S{z) and Si{z) in a systematic way beyound perturbation theory. The 



two kernels are defined in (42) and (p6|). Except for the first boundary vertex, 
they are formally the same. Therefore, w.l.o.g. we discuss in the following the 
kernel I^{z). Furthermore, as pointed out after Eq. we can set tg = 0. 



s < / — * 1 • ■ ^ • • — ► ^ < f 




Fig. 4. Illustration of the RG-method. Successively, the shortest contraction line is 
integrated out. The thick lines and dots indicate renormalized propagators and vertices. 



In section ^ we have found an effective theory in terms of the dot degrees of 
freedom. The reservoirs enter via pair contractions which couple the tunneling 
vertices. The aim is to integrate out all contractions in such a way that they can 
be interpreted as renormalization of Lq, G, A and B. The procedure is shown 
schematically in Fig. ^. We start with the integration over the shortest contrac- 
tion. This gives rise to a renormalization of G in Fig. ^. The integration over the 
next shortest contraction renormalizes the propagator e-'^^o{ti-t2) ^^j^g ^^q^ 
Proceeding in the same way, we find a renormalization of A and B in the next 
two steps. Finally we are left with a diagram which we can calculate easily by 
perturbation theory but with renormalized quantities. It can also happen that 
two or more vertices fall into one contraction which is integrated out, see e.g. 
Fig. ||. In this case we generate double-, triple-, and higher order vertices. In a 
perturbative renormalization group treatment one cuts this infinite hierarchy at 
a certain level. Here, we only consider propagator and single-vertex renormal- 
ization. 

As described, we first want to integrate over the short contraction lines, i.e. 
over short time scales of the function 7^(t). Short time scales correspond to large 
energy scales. This leads us to the usual way renormalization group is formulated. 
One introduces a cut-off function F{E/D) in the integrand of (pT|), where Z? is a 



An introduction to real-time renormalization group 15 




Fig. 5. Generation of a double-vertex. 



high-energy cutoff, F{0) = 1 and F{x) decays monotonically to zero sufhciently 
fast for — > ±oo. As a consequence, the pair contraction ^^^oit) depends on 
the cutoff D. In each renormaUzation group step one reduces the cutoff by an 
infinitesimal amount D D — 5D, i.e. one tries to integrate out a small energy 
shell. The mathematical challenge is to interpret the result of this integration as 
a renormalization of system parameters, like Lq, G, A and B. 

We prefer to develop the renormalization group procedure in real-time space. 
This turns out to be easier and more systematic. We introduce a cut-off function 
F{t/tc) which cuts off small time scales, i.e. F{x), with a; > 0, is a monotonically 
increasing function with F{0) — and lim^-^cx; ^"(a;) — 1, see Fig. ^. tc is a 
cutoff parameter for small time scales. We define a cutoff dependent reservoir 
contraction by 

l^^^tM ^ Ut) F{t/tc) ■ (75) 

If we choose a sharp cutoff function F{x) ~ 0{x — 1), ^^tA^) includes only those 
time scales which are precisely larger than tc- We note that the high-energy cutoff 
D introduced in Eq. (|2|) is independent of tc and is not used as a renormalization 
group flow parameter. Within our real-time formulation, it corresponds to the 
physical band width of the reservoirs. 

1 -- 

~ 1 / 



J >- 

X 

Fig. 6. The cutofT-function. 

Motivated by the picture presented at the beginning of this section, we hope 
that all time scales being smaller than tc, i.e. those being not present in j^t^{t) 
can be accounted for by renormalized quantities. To formulate this precisely let 
us write the kernel S{z) as a functional of Lq, G, A, B and 7 



S{z)=nL^,G%Al,Bl,^^{t)) , 



(76) 
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where p, /i,and t run over all possible values within the functional. After replacing 
7p by the cutoff-dependent function 7^4^ inside this functional, we try to find 
a cutoff dependence of all other quantities in such a way that the kernel stays 
invariant with the same functional J- 

S{z) - S,^ (z) + TiLot^ , G%,^ , AP,^ , B^,^ , 7^,^ (t)) . (77) 

This equation is only exact if we neglect higher order vertex corrections, as 
already pointed out above. Otherwise we have to include all higher-order vertex 
terms as arguments in the functional J-^ as well. Via 7^tt^, the second term on 
the r.h.s includes only time scales which are larger than ij,- AH other time scales 
are accounted for by the cutoff-dependence of Lg, G, A, and B. The first term 
on the r.h.s contains the contributions where all time scales are smaller than tc- 
This part is not included in the second term where at least one contraction line 
7/itc occurs. It is important to notice that St^{z) should not be viewed as being 



trivially defined by Eq. (77). If we choose a sharp cutoff function F{x) = 0(x— 1), 
St^ (z) can be contructively defined as containing only those diagrams where all 
contraction lines have a length smaller than tc- For tc 00 the second term on 
the r.h.s. of ( [77| ) vanishes, provided the function ^p,{t) decays sufficiently fast for 
t 00. Thus, we obtain the final solution from 

S{z) = lim St, (z) . (78) 

tc — ^00 

The aim is to find differential equations which describe the cutoff-dependence of 
IJtc{z), Lot,, Gtc, At, and Bt,. This will be the subject of the rest of this section. 

Let us increase the cutoff by an infinitesimal amount tc —>■ tc + dtc- Due to 
the invariance, Eq. (^ should not change 

= dE{z) + T{Lq + dLo, G + dG,A + dA,B + dB, 7 -I- ^7) , (79) 

where we have omitted the subindex tc and the indizes p and pi. The change of 
7 is known from the definition ( |75| ) 

dl^,tSt) = —i^{t)F'{t/tc)t/tl dtc , (80) 

or, for a sharp cutoff function F{x) — 0{x — 1) 

rf7p«cW = -l,.(tc)5{t - tc)dtc . (81) 



The increment d-y in (|79|) leads to all diagrams where one contraction line is 
replaced by ^7. We indicate this by a cross and define 



, d'y , 



G1...G2 = -{^Jl,{t,-t2)dtcGl{t^)...Gl,{t2), (82) 
We call this a "cross contraction". For a sharp cutoff function, we get 

G^^2 = 7^J (ic) S{t, -t2- tc) dtc Glit,) . . . G^:(t2) . (83) 
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Corresponding definitions hold for G replaced hy A or B. By convention we have 
included a minus sign into the definition since —d-y accounts for the time scales 
of the contraction between tc and tc + dtc, compare Eq. (pi]). We have to identify 
all terms created by d'y with a contribution arising either from dXJ, dLo, dG, dA 
or dB in Eq. (|7^) in order to fulfil invariance. Which one has to be taken depends 
on the number of vertices between Gi and G2, and wether one or both of the 
vertices are boundary vertices. 

Let us start with the simplest case of two successive vertices which are not 
at the boundary 

G1GX3G4 . (84) 

Here, the two vertices Gi and G4 are contracted to any other vertices. Motivated 
by Fig. ^, it is tempting to take this diagram together with G1G4 and interprete 
the cross contraction as a renormalization of Lq in the following sense 

G1G4 + / dt2dt3 G1G2G3G4 



1>2>3>4 

? 

J_ ^iLoti Qpi ^-i{Lo+dLo){ti-t4,) QP4, ^-iLoti {85) 

Expanding the exponential to linear order in dLg, we see that this requires the 
identity 

dt2dhGiG2G3G4 = -i / dtGi{dLo){t)G4 . (86) 

1>2>3>4 ti>t>t4 

To obtain such a relation we see immediately a problem. The integrand of the 
l.h.s contains a two-time object G2G3, whereas the integrand of the r.h.s contains 
the single-time object (dLo)(t). We have to decide to which time-variable we 
identify t. Let us make an arbitrary choice: t = t^. Then it is obvious to try the 
ansatz 

-i(dLo)(t3) = / dt2G2G3 . (87) 

2>3 

This is well-defined since it does not involve the time variables ti and t^. However, 
let us try to insert this ansatz into the r.h.s of (86). We get two terms 

-i J dhG4{dLo){h)G4= J d<2di3 G1GJG3G4 

1>3>4 1>3>4 

2>3 

dt2dt3 G1G2G3G4 + / dt2dt3 GiG2G3G4(88) 

1>2>3>4 2>1>3>4 

The first term agrees with the l.h.s. of (^6|), whereas the second one is a correction 
term which has to be considered separately. Let us try to interpret it as a term 
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arising from dG. We define 

{dGl\)'^'-\h) ^ - J dhd^GiChG^, (89) 



2>1>3 



wliere the subindex (1) indicates that this is the first contribution to dG (another 
one will follow below). Using this definition we finally get 

dt2dhGiChG^Gi ^ ~% f dt3Gi{dLo){t3)G4 + {dGP\Y^\h)Gi 



1>2>3>4 1>3>4 

dt2dUGiChG:iGi . (90) 

2>1>4>3 

In the last term on the r.h.s. the two time- variables ti and are "clustered" 
together, they both have to lie within the time interval [t3,i2]- Therefore, we 
interpret this term as a double- vertex which we neglect. In conclusion, we have 
achieved our final goal: we have interpreted a term arising from d'^ by a renor- 
malization of La and G. 

Let us try to understand what we have done so far and how to find a sys- 
tematic procedure for obtaining the complete RG equations without going again 
into the details of the foregoing exercise. We have defined the renormalization 
of Lq in (^^, which is natural since this expression contains only operators of 
the dot. We write ^ 

-z(rfLo)(f2) = G1G2 , (91) 

. . . . 

where the arrow indicates the time variable which is not integrated out. It is 

this time variable which is used for the definition of the interaction picture of 
dLo and, most importantly, which is used for the definition of time-ordering, see 
the r.h.s of Eq. (^6|). The time- variable ti no longer appears explicitly because 
it is an internal integration variable within the definition of dL{). This has the 
consequence, that time variables of other vertices do not care about the value of 
ti regarding time ordering. E.g., if we multiply {dLQ){t^) with a vertex G2 from 
the left, time-ordering only requires ^2 > ^3 but the ordering of t2 with respect 
to the internal integration variable ti is not prescribed. Thus, due to dLo, the 
following term will occur on the r.h.s. of ( [79| ) 

G2G^3 , (92) 

T T 

where ti > t2 > t^. However, this term is not present in the original series 
because the ordering of operators does not agree with the ordering of time vari- 
ables. Consequently, we have to subtract this term. This is the basic reason for 
the occurence of the second term on the r.h.s. of (|9C|). 

We interpret minus the correction term ( p2[ ) as renormalization of G. Again 
we have to specify the time variable for the interaction picture and for time- 



ordering. We choose t2 in analogy to Eq. (89). Together with the "real" vertex 
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renormalization, arising from a free vertex inside a cross contraction, we get the 
complete renormalization of G 

{dG%){t2) = G1G2G3 - G2G^3 . (93) 

We can now proceed to look at other correction terms which can occur by mul- 
tiplying (|9l| ) with two vertices from the left or (^) with one vertex from the left 
or right 

r4<n I ¥i — I rXn 1 X — 1 rXn 

G2G3G1G4 + G2G1G3G4 — G2G3G1G4 + G1G2G4G3 — G2G1G4G3 , (94) 

with ti > t2 > > ti. We have indicated the time-ordering variables of all sub- 
clusters which lead to these expressions. We see that the first and third term 
cancel each other. The other three correction terms have to be subtracted again. 
Minus the fifth term corresponds to the third term on the r.h.s. of Eq. (pO|). 
The second and fourth term are correction terms arising from the "real" ver- 
tex renormalization. However, all these correction terms correspond to double- 
vertices which we neglect consistently. Another double-vertex term arises from 
two free vertices within a cross contraction. Proceeding further in the same way, 
triple- and higher order vertex terms will occur which again are not considered 
here. We see that the procedure is very systematic and straightforward. 

In the same way we can proceed for boundary vertices. Since we want to 
calculate the Laplace transform S{z) of the kernel each diagram of the kernel 
(^2|) gets an additonal factor e*^*. Consequently we define the interaction picture 



for the boundary vertices slightly different from (41). We simply include the 
factor from the Laplace transformation 

Al{t) = e^'Ale-'^"' , BP{t) = e'^°*BPe~"' . (95) 

Analog to the above procedure we can write down immediately all terms with 
cross contractions containing boundary operators. They can be interpreted as 
renormalization of IJ, A and B 

dS{z) = iS2|t.=o , (96) 

{dAl\){t2) = A1G2G3 - A2ChG3 , (97) 

T T 



{dBP-){t2) = G1G2B3 . (98) 
T ^ 

Again, we integrate implicitly over all time variables which are not indicated by 

I X I I )^ I 

an arrow, with ti > t2 > t^. Terms like A1G2 or G1B2 do not occur since they 
do not lead to connected diagrams. For this reason, there is also no correction 
term in the renormalization group equation for B. 

Eqs. (|9l|), (p3|), and (|96|)-(|9^) are the final RG-equations. They can be trans- 
lated easily to get explicit expressions. In order to get the renormalization of the 
bare quantities we have to set ^2 = in all equations. To account for possible 
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(a) (b) (c) 

Fig. 7. Diagrams where additional minus occur if the ends of the indicated contraction 
are connected. In (a) a minus signs occurs if the particle number difference Ns — N^i 
is odd. In (b) and (c) a free vertex is crossed by connecting the end-points without 
crossing over po. 



minus signs arising from commutation of fcrmionic reservoir field operators, we 
have to include two auxiliary sign functions. They arise because we have to con- 
nect the two end points of the cross contraction. If the cross contraction couples 
the upper with the lower propagator, a minus sign occurs if the intermediate 
double-propagator state at the vertex with the larger time is odd (i.e. if the 
difference of the fermionic particle numbers on the upper and lower propagator 
is odd), see Fig. 0a. We account this by a sign operator af^ multiplying each 
cross contraction from the left. The matrix elements of this operator are defined 

by 

' (pp' for N, - N,, = odd 

[<y^ )ss'.,s' - I ^ for iV, - iVy = even ' 

where TV, is the fermionic particle number of state s and /it is a fermionic con- 
traction (for bosonic contractions no additional sign has to be considered). The 
second sign function concerns the case when one free vertex occurs within a cross 
contraction. A minus sign occurs when the cross contraction crosses over the free 
vertex when we want to connect the two end points. Since a is multiplied from 
the left by definition, we connect the vertices always in such a way that we do 
not cross over po. We denote by p, pi, and p[ the indices of the free vertex, the 
vertex with the larger and the one with the smaller time variable of the cross 
contraction, respectively. If p = p[ we get an additional minus sign, see Fig. |^ 
and Therefore we multiply the "real" vertex correction by an additional sign 
function 

pp'i ^ f -PP'i for ^J■ and ^' fermionic 

'mm'i I 1 otherwise ' ^ ' 
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oo J 1 





oo , 
— oo 










oo , 
— oo 






; Ty^A a^/'^ AP^X mi - APy,\''^ G^X (t) ' 


G^,{t') 


oo / / / / 
-oo 



(101) 



Implicitly we sum over all double indices on the r.h.s. of these equations which do 
not occur on the l.h.s.. We note that the overall sign on the r.h.s. differs from the 
one in [Q since we have included the factors —i into the vertices. Furthermore, we 
have generalized the RG equations to an arbitrary cutoff- function F{t/tc). For a 
specific cutoff-function, e.g. F{x) — 0{x — 1), and taking matrix elements of the 
RG-equations, all integrals in ( |101| ) can be calculated analytically. We are left 
with pure differential equations which can be solved numerically in a straightfor- 
ward and very efficient way. Finall y, t he asymptotic value limt^_>oo ^tS^) gives 
the solution for the kernel, see Eq. (|7q). Using (46) and (47), we get the complete 
time-evolution of the reduced density matrix of the dot. The initial condition 
for the boundary operators at ic = is given by = B'p^ = G^. To get the 
current kernel Z'/(z), we take the same RG-equations with the initial conditions 



= II and Bl 



Gl- 



Numerically, one can not take the initial ffow parameter t^^ — Q and the 
final one tl oo. However, one can check that the final solution is stabel for 
sufficiently small i° and large . This means that l/tj? and l/t^ have to be much 
larger resp. smaller than all other energy scales 



^ < Tl, Tfl, T, eV, |e| < D < ^ 



t: 



(102) 



where eV = fiL ^ I^r is the bias voltage, and Fi^jFj^ are defined in (|23|). In 
addition, we want to send the band- width cutoff D ^ oo. This means that we 
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have to check the stabihty of the sohition for ^ 0, D ^ oo, with Dt^ 0. 
For the problem under consideration this is indeed possible. 

There are three essential differences of our final RG equations to conventional 
poor man scaling and operator product expansion techniques [ p"2[jl^ : 



1. We have formulated the RG-equations within a real-time formalism 
on the Keldysh-contour. 

2. We have not expanded the interaction picture of the vertex 
operators for small times, i.e. the renormalized propagation exp(±zLoi) 
is taken fully into account. 

3. For tc > t° we get Lq ^ [Hq,-], i.e. we generate non-Hamiltonian 
dynamics during RG. 



The first point is necessary for describing nonequilibrium phenomena. The sec- 
ond one is important for self-consistency reasons (Lq in the exponent is a renor- 
malized quantity) and for the stability of the solution for t[ ^ oo. Only if one 
expands the exponentials exp(±zLo^) in ti one has to stop the RG-flow when 
\Xtc\ ~ 1, where A denotes any eigenvalue of Lq. We do not need this expansion 
and, consequently, are able to integrate out all time scales. The third property 
is the most important one. It is essential for a nonequilibrium theory to describe 
the physics of dissipation. Thus, the renormalized supcropcrator Lq should no 
longer be expressible by a commutator with a renormalized Hamiltonian _ffo, i-e. 
in matrix notation 

iLo)siS2,s[s'^ ^ iHo)siS2'>s{s'2 - SsiS2iHo)s'2s{ ■ (103) 

To show this let us first define the renormalized Hamiltonian in a natural way. It 
is obvious that contractions which connect the upper with the lower propagator 
of the Keldysh contour, do not lead to Hamiltonian dynamics. This was already 



explained at the end of section 2.2. Thus, in order to define the renormalized Hq 
let us consider the same RG equations but allow only for contractions within the 
upper or lower propagator. Formally, this means p = p' and pi — p'l in ( |101| ). 
Under these conditions we can write the solution as 

{L0)S^S2,S[4 = {Ho)s,S2Ss'^s'^ - Ss,S2{HI)s'^s'^ ■ (104) 

To obtain this, we have used the symmetry relations (|5^) together with the 
property 

= (105) 

The renormalized Hamiltonian Hq is defined by considering only the upper prop- 
tor, i.e. settings = p' = pi = p'^ = + and replacing Lq — > Hq and G+ —ig^i 
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in the RG-equations. We note that Hq is non-hcrmitian since we generate com- 
plex energies during RG. Physically this describes broadening of levels or fi- 



nite life-times. From (104) we conclude that the backward propagator involves 



Hq , i.e. the sign of energy-broadening is different for the forward and backward 
propagator. This means, that even if we disregard contractions connecting the 
upper with the lower propagator, we get Lq 7^ [-ffo,"]- However, this does not 
imply the physics of dissipation. It describes the physics of a finite broadening 
of the energy levels due to the coupling to the environment. This reflects the 
Heisenberg-uncertainty relationship. 

In contrast, if we couple the forward and backward propagator by contrac- 
tions, the matrix Lq will completely change its form. It can neither be represented 



as ( |103D nor as (104). The couphng of the propagators gives rise to the generation 
of rates where the states on the upper and lower propagator are changed simul- 
taneously. Rates describe the evolution of the system into a stationary state, i.e. 
we generate irreversibility or dissipation during RG. 

The RG equations preserve conservation of probability and give real expec- 
tation values for hermitian observables. Once the symmetry relations and sum 
rules, stated in Eqs. (^-(p4|) and (|60|)-(|6l|), are fulfilled initially, they are not 



we 



changed by the RG-flow. Thus, by applying the same proof as in section 3.1 
obtain a normalized probability distribution of the dot and a real expectation 
value for the current. 

Under certain circumstances the RG equations can be simplified. This will 
be important for the exact solution presented in the next section. We decompose 
the pair contraction (^0|) trivially into two terms 

7p(0-7'(i)+7p(0 , (106) 

with 

lt^l{[j^.{t),J^U) , % = \{[jf.it),J,]a) , (107) 

where cr = ± for /i = bosonic (fermionic), and [•, ■]^a- denotes the (anti-)commu- 
tator for cr = ±. This decomposition is useful since in many problems, the part 
7^(t) turns out to be proportional to a delta function or derivatives of a delta 
function (at least for the band width cutoff D oo). For our special model we 
get (compare (p6|)) 

7^(0 = Y<^(t) ■ (108) 
Such a contribution can be included into the initial conditions of the RG equa- 

into (0) 

Dr>'/,\ r / - /,\n f 1 for u = bosonic 

iZ (i) - [HW + \ y for = fcrnriomc ' ^''^^ 



tions. First, we decompose 7|^^/(0 = 7^^'(i) +7^^'(0 ^ '^sU by inserting (|10 



We have obtained the explicit dependence on the indices p and p' in this equa- 
tion. Therefore we try to include these factors into the definition of the vertex 
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operators. Having included the 7^-parts into the initial conditions, it turns out 
that closed RG equations can be found for the quantities 

_ ( J2piGfi)ss',-- for Ns ~ Ns' = even or fi = bosonic 



^ 'iZ^pPCCfJss', - Otherwise 

(110) 

with Ns being the fermionic particle number of state s. In the same way we 
define the boundary vertices and B^. We have included a factor —i for the 



second case in (110) in order to get the same symmetry relations as in ( |5l| ) 

(Ga')siS;,S2S2 ^ (^p)s'iSl,S^S2 ' (Hi) 

Using the fact that fermionic vertex operators {Gfj_)ss' .sis[ change the parity of 
the fermionic particle number difference A^s^ — Ng'^ Ng — Ns' , we get after a 
lengthy but straightforward calculation the following RG-cquations arising from 
the 7^-parts 







00 

^G^ = -Jdt J dt'!^{t~t')[a^^,G-,At)G, - G^Gf^At)] G^.(t') 

-00 



00 

^^A^^-Jdt J dt'^{t~t')[a^^,Af,,{t)G^ ~ A^G-^,{t)] G^At') 
-00 



cx; 



B^ = -Jdt J dt'q^it-f)a^^,Gf,,it)GMt') 



(112) 



where 



and 



~ - /.N J 1 for H = bosonic , , 

7,(t) -7,(i)^ i^,^onic ' (1^3^ 

^ f 1 for /i or = bosonic (114) 

I ~1 I'-'r M ^-nd fi' = fermionic 

We see that the new RG equations are more compact but we note that they 
can not be applied to all models. Sometimes, like e.g. in spin boson models, it 
happens that the part 7^(t) can not be expressed by the derivative of a delta 
function and, consequently, can not be incorporated into the initial conditions. 
In this case, one should use the original and generally valid RG-equations (|lOl|). 
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5 Exact solution of the RG equations 



In this section we demonstate that the RG equations can be solved exactly for 
the special model (|l|)-(||) of a single non-degenerate dot state. It turns out that 



the result for the kernels agrees with that of section 3.2. We conclude that the 
RG-approach solves the present model exactly. 

We ch oose the sharp cutoff-function F{x) — 9{x — 1) and get from the RG- 
equation (101) for a matrix element of the kernel 



dUlfAt,)e-'^ (4), 



[Bl )s.-s.rs-s • (115) 



Here, due to particle number conservation, we have used the fact that the only 
nonvanishing matrix elements of Lq are as = (io)ss,ss and {Lq)ss,s' s' ■ To eval- 
uate this equation we need the tc dependence of as and the boundary vertex 
operators. To obtain this it is convenient to include the part 7^ into the ini- 
tial conditions. The latter follow from inserting the part 7^^ {tcY into the RG- 
equations (101), and integrating. We denote the initial value for Lq by Lq. The 



vertex operators do not obtain any initial renormalization from the (5-function 
parts since there is no phase space for the time integrals in (101). Taking the 
part {tcY from the first term of ( |109| ) together with (108), and inserting it 
in the RG-equation (101) for Lqi we find 



^^5{tc)pp'{GlGl)s-s,s-s 



This gives 



r 



r 



Qfi = e — I- 



(116) 



(117) 



- 2 ' ^ 2 ' 

where the part with the single particle energy e stems from the initial value 
without any renormalization. 

As we will show in the following we find no further renormalization from the 
7^-part, i.e. using the simplified set of RG-equations (112), we get 



dt, 



-as 



dt, 



dt 



(118) 



This means that the renormalization of the kernel due to the 7^-part can be 
calculated with unrenormalized vertex operators and using the result (117) for 
as- Since, trivially, the same applies to the 7^-part, we can directly evaluate 
(115) in this way and get the result ( 7^ of section 3.2. 

What remains to be shown is Eq. (118). We first note the following properties 



of the matrix elements of the vertex operator, which follow from the definition 

^ ^j jGr fj,) ss^s' s' — , {G ss s' — {G ss^s' s' ■ (ll'^) 



(110) 
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Furthermore, using (|53|), we get 



— iLnt 



= 1 . 



(120) 



We apply these properties to the foUowing expression 

/ , V^^ IJ. ) ss.sisi ye jsisi,S2S2 V"-^P 

SlS2 

*e~*" = '* (6*^)53,5151 ^(Gp')s2S2,s'5' =0. (121) 



-■p' JS2S2,S'S' 



Using this result after taking matrix elements of the RG-equations (112), we find 
immediately (|ll^). 



We conclude that after having set up the general RG-equations (101) and 



(112), the exact solution of the present model can be found analytically in a 
straightforward way. There is no need to consider diagrammatic details as in 
section 3.2. Furthermore, we see that the effect of the reservoirs is simply a 



broadening of the local state of the dot by r/2, see (117). Otherwise the evalu- 
ation of the kernel is the same as in lowest order perturbation theory in F. 



6 Summary and outlook 

In this article we have presented a new viewpoint to analyse nonperturbative 
aspects of nonequilibrium systems. We considered a small system coupled lin- 
early to several baths. Nonperturbative means that the coupling between system 
and bath is so strong that quantum fluctuations induce broadening of the states 
in the system together with possible renormalization of energy levels and the 
coupling to the environment. In macroscopic systems, such effects are negligible 
since the interaction with the environment is a surface effect. In contrast, our 
goal is to describe mesoscopic systems like quantum dots, magnetic nanoparti- 
cles or chemical molecules coupled to particle or heat reservoirs. Usually such 
systems are treated within perturbation theory ||lj,|l5 16 . This means that the 



kernel S(t) of the kinetic equation is calculated in lowest order perturbation 
theory in the coupling, the so-called golden rule or Pauli-Master equation ap- 
proach. Our aim here was to find a systematic way to consider self-consistently 
an infinite series of higher-order contributions to the kernel. This is a nontriv- 
ial task, especially in nonequilibrium where a real-time formalism on a Keldysh 
contour has to be used. Renormalization effects known from equilibrium theories 
have to be incorporated consistently within a kinetic equation. Our point of view 
relies on renormalization group ideas. We try to integrate out all energy scales 
of the bath in infinitesimal steps. Each step is interpreted as a change of various 
quantities, like broadening and renormalization of energy levels, and generation 
of rates. For formal reasons we have set up this procedure in real-time space. 
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During this procedure we keep the kernel of the kinetic equation invariant and 
generate non-Hamiltonian dynamics to describe the physics of dissipation. The 



final RG-equations are presented in ( |101| ), and, alternatively, in (ff2). Solving 
them, provides the complete description of the time-evolution of the reduced 
density matrix of the system and arbitrary observables being linear in the field 
operators of the bath. This includes the consideration of an initally out of equi- 
librium state as well as the description of stationary nonequilibrium situations. 
We have solved the equations exactly for the special case of a quantum dot 
with one state coupled to two particle reservoirs. The solution turned out to be 
identical to the exact one. 

The reader might argue that this is a trivial result since the model under 
consideration is very simple. However, as shown in a recent paper, the same 
RG equations provide not only the exact solution for a noninteracting quan- 
tum dot with one state, but gives a good solution also for the spin-degenerate 
case including a finite on-site Coulomb interaction U , the so-called Anderson 
impurity model in nonequilibrium. Here, not only energy broadening but also 
energy and coupling constant renormalizations are important. In the mixed- 
valence and empty-orbital regime, i.e. for level positions near or above the Fermi 
level of the reservoirs, it was shown that the linear conductance and the average 
occupation in equilibrium agrees perfectly with Friedel sum rules, Bethe ansatz, 
and numerical renormalization group methods within 2-3%. This means that 
the RG provides a good solution for the whole range from U = to ?7 = cx). 
This is a surprising result since usually methods designed for strong interaction 
do not work well for weak interaction and vice versa. Therefore, the fact that 
the RG gives the exact result for the noninteracting case is not at all a trivial 
result. In fact, within the slave boson technique [ p^ , which is a well-known and 
well-established method in the theory of strongly correlated Fermi systems, it is 
very complicated to find the exact solution for U — . 

Another example where the RG-method has been apphed is the study of 
transport through the metallic single-electron transistor ||l8[. For this case, the 
RG-equations have been solved in sixth-order perturbation theory in the cou- 
pling. Surprisingly, it turned out that the solution agrees very well with exact 
perturbation theory in the same order. The reason why the RG-equations have 
not been solved in all orders is that the result was not finite for the band width 
cutoff D — > oo. Only if certain parts of double-vertex corrections were included, 
the solution turned out to be cutoff-independent. A more detailed study of the 
D-dependence and the solution for higher orders is currently under way. 

The RG-equations on the level of propagator and single-vertex renormaliza- 
tion are pure differential equations and, consequently, can be solved numerically 
very quickly. Typical response times range from seconds to a few minutes for one 
set of parameters, even on a usual PC maschine. There are still some problems 
with the asymptotic solution of the RG equations for tc ^ oo. First, it is not 
known rigorously wether a stationary solution exists at all. Secondly, the numer- 
ical solution is plagued by oscillating functions, typical for real-time problems. 
However, for the problems studied so far, and provided the coupling is not too 
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strong, there is a stationary numerical solution over a sufficiently long period in 
tc- This applies at least to the physical quantities under consideration, like the 
probability distribution and the current. 

Another reason for the efficiency of the present method is the possibility to 
calculate physical observables directly without the need of correlation functions 
like in linear response theory. Furthermore, if desired, correlation functions can 
also be studied with the RG-method. In this case, there are additional RG- 
equations to describe the renormalizations of the external vertices defining the 
correlation function. These equations together with explicit solutions will be 
published in forthcoming works. 

The method is also applicable to the study of the ground state energy since 
the RG-flow on the single forward propagator provides the S-matrix. This idea 
has been appHed to the single-electron bo x |l0[| , coupled metaUic islands 
and the one-dimensional Polaron problem [^0[. In the first two cases very good 
results have been obtained, even comparable to very time-consuming QMC- 
simulations. For the Id-Polaron problem, the results were at least satisfactory 
for small coupling. Here, problems occured since the correlation function of the 
bath does not decay for long times and undamped modes occur which make the 
numerical analysis of the asymptotic solution very difficult. 

Finally, we remark that a challenge for future research is the consideration of 
higher-order vertex corrections. Our RG-scheme provides a systematic treatment 
for setting up RG-equations for all kinds of multiple vertices. However, even on 
the level of double-vertices, the number of terms increases considerably and the 
RG-equations become integral-differential equations. The reason is the retarded 
nature of double- vertices. Therefore it is necessary to find physical arguments 
to select the most important terms or to improve the numerical efficiency by 
neglecting the retardation in a convenient way. Whereas single-vertices describe 
basically charge fluctuations (in case of coupling to particle reservoirs), double- 
vertices describe physical processes via virtual intermediate states where the 
local system can change its state without changing its particle number. This is 
e.g. important for the study of spin fluctuations in local impurities or quantum 
dots, leading to the Kondo effect ||^,|l^. Although such processes are pertur- 
batively included in the RG-equations set up in this article, it is important to 
study them fully sclf-consistently by considering corresponding RG-equations for 
double- vertex terms on a Keldysh contour. 
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Appendix 

Here we prove Eqs.(62)-(|4|). The first one follows from 

{Gf)ss.ss + {G^,,.,s - , (122) 



which can easily be seen by using the definitions (g7|) and (^8|). The reservoir 
contraction is the same in both terms since is the vertex at later time, 
compare The same proof can be used for 1^ by using the definition (|37|). 
Since there is no change of sign of 1^ if p is changed to p, we have to add a factor 
p under the sum in (^). 

To show the second property (^), we again use the definitions ( p7| ) and (p8|), 
and find | p= | p= 

(G^G^O.s-.. + (G^Gf )..-,■■ =0. (123) 

Both reservoir contraction are the same in both terms since both vertices are at 
later time. However, the reader can convince himself very easily that there is an 



additional relative sign between the two terms of (123) due to the interchange 
of reservoir Fermi operators. 

The third property ( |64[ ) follows from 

hm {(GpJ).,,. + (Gpf ),,,.} = . (124) 

Here, by comparing the two terms, the contraction associated with the second 
vertex changes from 7^ to 7^''. Otherwise the two terms are the same. However, 
the sum 77(i) + Jr^it) gives a ^(t)-function in the hmit D ^ 00, see (p^). This 
means that there is no phase space in time for the first vertex and the sum is 
zero. 
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